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We construct a family of time-independent nearest-neighbor Hamiltonians coupling eight-state 
systems on a ID ring that enables universal quantum computation. Hamiltonians in this family 
can achieve universality either by driving a continuous-time quantum walk or by terminating an 
adiabatic algorithm. In either case, the universality property can be understood as arising from an 
efficient simulation of a programmable quantum circuit. Using gadget perturbation theory, one can 
demonstrate the same kind of universality for related Hamiltonian families acting on qubits in 2D. 
Our results demonstrate that simulating ID chains of spin-7/2 particles is BQP-hard, and indeed 
BQP-complete because the outputs of decision problems can be encoded in the outputs of such 
simulations. 



I. INTRODUCTION 

With quantum circuits [I] 12 [3], one can decompose 
even the most complex quantum computation into a 
sequence of simple operations called gates that act on 
simple parcels of information called qubits. Mathemat- 
ically, the action of a T-gate quantum circuit on an 
rt-qubit pure state can be expressed as Up ■ ■ ■ Ui\yj), 
where l^) represents the input state and each Ui rep- 
resents a unitary transformation drawn from a quantum 
gate basis, such as the universal gate bases described in 
[H HI El El El [H El [IHl [HI H21 H31 HU • This mathematical 
description is usually interpreted as a set of instructions 
for applying the gates U\ , U2 , . . . , Up in sequence to the 
qubits. Feynman noted in Ref. |15j that this description 
can instead be interpreted as a blueprint for a device in 
which the gate sequence is "printed" directly into hard- 
ware. Specifically, from the mathematical description of 
a quantum circuit, Feynman constructed the following 
time-independent Hamiltonian that acts collectively on 
both the input state \yb) and an auxiliary system he called 
a "cursor" consisting of (T + 1) qubits: 
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where er^f and cr^~ denote raising and lowering operators 
on the fcth cursor qubit. Feynman proved that evolution 
for a time T/2 under Hp with the cursor initialized in the 
state |1, 0, . . . , 0) will efficiently approximate a sequenced 
implementation of the circuit. In the language of quan- 
tum information, one would say that Feynman proved 
that continuous-time quantum walks |16j are universal 
for quantum computation. 

Feynman's construction has received little attention as 
a possible quantum computing architecture. This is most 
likely because a direct implementation of Hp appears to 
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be far beyond the reach of any known quantum technol- 
ogy. Some of the difficulties are that Hp has a) four-body 
interactions (when two-qubit gates are used), b) spatially 
nonlocal interactions (because the cursor gets far away 
from the data as T and n grow), and is c) exponentially 
sensitive to decoherence via Anderson localization [17] 
(because the quantum walk generated by Hp is effec- 
tively on a uniform line). The central goal of this paper 
is to construct a Hamiltonian that achieves the same task 
as Feynman's but which is simpler and hopefully closer 
to a technologically-reachable implementation. 

There is a wealth of previous work on simplifying Feyn- 
man's Hamiltonian in the context of exploring the com- 
putational universality of adiabatic algorithms and of 
demonstrating the QMA-complcteness of finding Hamil- 
tonian ground states PH [E0 HH Effl H3 [23] . (Actually, 
these efforts are directed at simplifying a variant of Feyn- 
man's Hamiltonian proposed by Kitaev in |18j.) How- 
ever, these simplifications are done with the goal of find- 
ing Hamiltonians whose spectra are close to that of Hp 
but which are not necessarily close to Hp dynamically. 
This latter kind of closeness is achieved when the oper- 
ator norm of the difference between the Hamiltonians is 
close, because, as the following well-known inequality for 
Hermetian operators indicates [24], this implies that the 
corresponding dynamics generated by the Hamiltonians 
is close: 
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Examples of Hamiltonians whose spectra can efficiently 
be made close to that of Hp include a nearest-neighbor 
Hamiltonian on a two-dimensional grid of qubits |21j . 
a translationally-invariant nearest-neighbor Hamiltonian 
on a line of 30-level systems [22], and nearest-neighbor 
Hamiltonians on a line of 9-level systems 23J. 

Although consideration of Feynman's quantum com- 
puting architecture has been hampered by technologi- 
cal unreachability, it has not been overlooked entirely. 
For example, Margolus [55] presented an early gener- 
alization in which Hp is replaced by a spatially ho- 
mogeneous finite-range 8-body Hamiltonian on a two- 
dimensional torus of qubits. Spatial homogeneity makes 
this model particularly simple — it might best be termed a 
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continuous-time quantum cellular automaton [22 . 2"S1 [2T1 
28, 29 . However, the 8-body interactions and torus ge- 
ometry keep it from being any more practical than Feyn- 
man's original architecture. More significantly, Margo- 
lus' model is only universal for classical, not quantum 
computation. Janzing and Wocjan [30 improved Mar- 
golus' model so that it is programmable, universal for 
quantum computation, and homogeneous on the surface 
of a cylinder, but at the cost of requiring 10-body inter- 
actions among the qubits over a finite range. Janzing |31j 
has further modified this construction so that with mild 
spatial inhomogencitics in the interactions, it achieves 
programmable universality using only nearest-neighbor 
interactions among 3-level systems on a two-dimensional 
lattice. Recently, Nagaj and Wocjan removed these inho- 
mogeneities in work concurrent to ours |32j . and demon- 
strated that programmable universality can be achieved 
by a continuous-time cellular automaton in one dimen- 
sion; their construction requires increasing the number 
of levels per system from three to ten to achieve this. 



In this paper, we present a programmable, universal 
architecture for quantum computation that uses a time- 
independent nearest-neighbor Hamiltonian on a one- 
dimensional ring of 8-lcvel systems. When combined with 
the results of Oliveira and Terhal about the operator 
norm closeness of Hamiltonians generated by a certain 
"gadget" perturbation theory technique |21J . we show 
that we can modify our architecture to use only qubits at 
the expense of requiring a two-dimensional rather than 
a one-dimensional geometry. Additionally, we show how 
our model can also be used to reduce the state space per 
carrier from 9 to 8 in a recent proof of the universality of 
the adiabatic algorithm [23.. 



The remainder of our paper is organized as follows. 
In Sec. |III[ we motivate our architecture by defining a 
relevant programmably universal quantum circuit family 
T. In Sec. |IV| we define a family of one-dimensional 8- 
state Hamiltonians on a ring whose interactions mimic 
the gates used by circuits in T . In Sec. |IV| we show how 
these Hamiltonians can be used to drive continuous-time 
quantum walks that simulate circuits in T efficiently. In 
Sec. [V] we demonstrate the universality of our architec- 
ture under adiabatic evolution, following similar argu- 
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ments to those of Aharonov et al. [20] [23] . In Sec 
present a method for realizing this model using nearest 
neighbor Hamiltonian interactions between qubits on a 
regular rectilinear lattice on the surface of a cylinder by 
leveraging results of Oliveira and Terhal [21]. Sec. VII 
concludes, followed by Appendix [A] which steps through 
how an example 3-gate simulation would be effected in 
our model. 



II. A PROGRAMMABLY UNIVERSAL 
QUANTUM CIRCUIT FAMILY 

A. Notation and basic definitions 

We begin by reviewing relevant notation and basic def- 
initions, much of which can be found in standard text- 
books on quantum computation, such as [18, 33 . A qubit 
is a quantum two-state system, the computational basis 
of which we define as |0) := (J), |1) :— (°). A gate is a 
unitary transformation on qubits and a circuit is a com- 
position of gates. Gates discussed in this paper include 
the NOT gate, the Hadamard gate, the phase gate, and 
the swap gate, defined in the computational basis as 
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Each gate may be depicted as a circuit element in 
which qubits are represented by horizontal lines (wires) 
entering on the left and exiting on the right. The circuit 
elements depicting the gates above are 
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The controlled-U gate, denoted A(t7), applies the gate 
U to a target state conditioned on a control qubit being 
in the state |1). The matrix and circuit element repre- 
senting A(U) are 



A(U) := 



/ 
u 
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where / represents the identity matrix with dim / = 
dim U . 

We use two non-standard shorthand notations for gates 
in this paper. The first is that we denote gates controlled 
by the state |0) rather than the state |1) as 



U 



U 



The second is that we denote a cascaded array of nearest- 
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neighbor swap gates as 



A gate basis is a set of gates; a circuit is expressed over 
a gate basis Q if all of its gates are elements of Q. A cir- 
cuit family is uniformly generated if a description of its 
elements can be constructed by finite means, e.g., by a 
Turing machine [34 . A gate basis Q is universal if for any 
unitary transformation U and for any desired precision e, 
there is a uniformly generated circuit family over Q that 
can approximate U to within e. An example of a uni- 
versal gate basis is the Kitaev gate basis {H, A(5)} [13] . 
A programmable circuit is one that can be conceptually 
divided into "program" and "data" regions, where the 
state of the program region controls which gates are to 
be applied to the state of the data region. For example, 
Fig. [I] depicts a programmable circuit over the gate ba- 
sis {A(H), A 2 (S)} in which the program controls which 
gates from the Kitaev basis are to be applied to the data. 
A gate basis Q is programmably universal if for any de- 
sired precision e and any number of qubits n, there is a 
uniformly generated programmable circuit V over Q such 
that for any n-qubit unitary transformation U , there is 
a uniformly generated program state that enables V to 
approximate U to within e on the data portion of V. 
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FIG. 1: A simple programmable circuit over the Kitaev gate 
basis. 



B. Programmably simulating circuits over the 
Kitaev gate basis exactly 

Consider the programmable circuit V depicted in 
Fig. [2] Its program region is decomposed into three parts 
labeled cursor, swap, and Hadamard. Given a description 
of a quantum circuit U on n qubits containg T gates that 
is expressed over the Kitaev gate basis, the program re- 
gion of V can be initialized so that 0(nT) iterations of 
V will simulate U exactly. To see how, we first trace 
through one pass of the circuit V to see how it responds 
to various initializations of the program regions. We then 
describe how to initialize these regions to carry out the 
simulation of U. 
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FIG. 2: The quantum circuit V, whose action is described in 
the text. 



Phase I. The £ gate moves the upper-most swap qubit 
state to the lower-most swap qubit and moves every other 
swap qubit state up one qubit to accommodate. 

Phase II. If the lower-most swap qubit is in the state 
|0), then the two-qubit cursor state permutes from 
j 00) -> 1 00) and from |01) -> |10) -> |11) -> |01). 
We depict this permutation in glyphs as (3) — > (3) and 
(J) — > © — > O "~ * ®- The reason for this choice of 
glyphs will become apparent later. 

Phase III. The £ gate moves the upper-most data qubit 
state to the penultimate data qubit and moves the state 
of every data qubit in between up one to accommodate. 

Phase IV. If the two-qubit cursor is in the state (j), 
then the penultimate data qubit state is swapped with 
the final data qubit state. If the cursor is in the state 
O, then first A(5) is applied between the last two data 
qubits, followed by an application of H to the last data 
qubit, conditioned on the upper-most Hadamard qubit 
being in the state |1). If the cursor is in any other state, 
no action is performed on the data qubits. 

Phase V. The S gate moves the upper- most Hadamard 
qubit state to the lower- most Hadamard qubit and moves 
every other Hadamard qubit state up one to accommo- 
date. 

Given this action of one iteration of V, we program V 
so that it follows a sequence of "behaviors" that depend 
primarily on the state of the two-qubit cursor: 

Behavior 1. The two-qubit cursor is initialized in the 
state (j) and the first n — k swap qubits are initialized in 
the state |1). The effect of applying V a total of n — k 
times given this initialization is to cyclically permute the 
states of qubits in each program region and in the data 
region upwards n — k times. Importantly, the state of 
the fcth data qubit from the top before this behavior is 
stored in the lower-most data qubit after this behavior. 
The glyph (j) is intended to remind that when the cursor 
is in this state, the data qubit states are all cyclically 
permuted. 

Behavior 2. The cursor is initialized in the state ®, 
the upper-most swap qubit state is initialized in the state 
|0) and the next n — j — 2 swap qubits are initialized in 
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the state The effect of applying V once given this 
initialization is to cyclically permute the swap qubit's 
states upwards once, cycle the cursor state from (J) to 
®, cyclically permute all but the state of the last data 
qubit upwards once, and cyclically permute all Hadamard 
qubit states upwards once. The effect of applying V a to- 
tal of n — j — 2 more times given this initialization is to 
repeat this behavior n — j — 2 times, except it will not cy- 
cle the state of the cursor qubits. Importantly, the state 
of the jth data qubit from the top before this behav- 
ior becomes stored in the third-to-last data qubit after 
this behavior. The glyph @ is intended to remind that 
when the cursor is in this state, the effect is to cyclically 
permute all but the last data qubit state. 

Behavior 3. The cursor is initialized in the state ® and 
the upper- most swap qubit is initialized in the state |0). 
The effect of applying V once given this initialization is 
to cyclically permute the swap qubit states upwards once, 
cycle the cursor state from @ to Q, cyclically permute 
all but the last data qubit state upwards once, apply 
A(S) between the last two data qubits, apply H to the 
last data qubit if the upper-most Hadamard qubit is in 
the state 1 1) , and cyclically permute the Hadamard qubit 
states upwards once. Importantly, the states of qubits k 
and j manipulated in behaviors 1 and 2 have the gate 
A(S) applied between them and the state of qubit k has 
the gate H applied to it conditioned on the state of the 
upper-most Hadamard qubit before this behavior. The 
glyph O is intended to remind that when the cursor is 
in this state, the effect is to apply gates from the Kitaev 
gate basis to the data qubits. 

Behavior 4. The cursor is initialized in the state O an d 
the upper-most swap qubit is initialized to the state |0). 
The effect of applying V once given this initialization is to 
cyclically permute the states of the swap qubits upwards 
once, cycle the cursor state from Q to (J), and cyclically 
permute the data and Hadamard qubits upwards once. 
Importantly, after this behavior, the cursor is in the same 
state it started in for behavior 1. 

Behaviors 1-3 describe how, given a suitable initializa- 
tion of program qubits, the gate (I<g)H)A(S) or A(S) can 
be applied to any two desired data qubits. Moreover, be- 
havior 4 shows that by lengthening the program region, 
any sequence of such gates can be implemented. How- 
ever, the Kitaev gate basis is not {A(5), (J <g> H)A(S)}, 
but {H, A(5*)}. In order to apply H to any desired qubit, 
we use a trick: we program V to apply (7 ® H)A(S) ■ 
A(S) 3 . Because A(S) 4 — I, this is the same as applying 
H. 

If a circuit U over the Kitaev gate basis contains Ia(s) 
controlled- S gates and Th Hadamard gates on n data 
qubits, then the number of program region qubits (and 
gates) needed by V is at most 4n(TA(s) + 2Th), so the 
simulation by V is linear in the size of U. 

It is worth noting that V is defined over a gate ba- 
sis that contains gates that are not spatially local when 
the qubits are arranged in a line. Spatial locality in one 



dimension could have been achieved much more easily 
(and without the need for programmability) with, say, 
the gate basis {H, A(S), SWAP}. Nevertheless, all the 
nonlocality in V is between the cursor program qubits 
and the other data qubits and enables V to be repeated 
in a simple fashion to simulate arbitrary circuits U over 
the Kitaev gate basis exactly. By replacing each of the 
non-cursor qubits with eight-state systems (to simulate 
one data qubit plus two cursor qubits), it is possible to 
derive a related programmable circuit that requires only 
spatially local gates in one dimension. Instead of going 
through this exercise, we move on to how to construct 
a Hamiltonian that would simulate such a circuit via a 
continuous-time quantum walk. 



III. A FAMILY OF NEAREST-NEIGHBOR 
HAMILTONIANS COUPLING 8-LEVEL 
SYSTEMS ON A RING 

In this section, we construct a Hamiltonian H s acting 
on nearest-neighbor eight-state systems on a ring. The 
interactions in H s are closely related to the gates used by 
the circuit V of Fig. [2] of the last section. We defer a dis- 
cussion of how Hg, is used to achieve universal quantum 
computation to subsequent sections. 



A. State space and geometry 

Let V be the circuit depicted in Fig. [2] and let V^ K \ 
K = G{nT), be a circuit with suitably initialized inputs 
that simulates a T-gate, n-qubit circuit over the Kitaev 
gate basis. Let 7ig denote the Hilbert space of a quan- 
tum system having eight possible orthogonal states and 
let Tif K denote the Hilbert space of a collection of K of 
these systems arranged in a one-dimensional ring. For 
convenience, we will think of the eight-state system as 
arising from the combined state space of quantum sys- 
tems having two and four states respectively. In other 
words, we will use the conceptual decomposition of the 
Hilbert space of each of these systems as "H% = 'Hi®'H^ 
where Tld indicates a Hilbert space of dimension d. 

As illustrated in Figure [3] this ring can be divided 
into two parallel "lines" and three regions. This layout 
mimics that of the circuit V rotated on its side, but with 
periodic boundary conditions. However instead of having 
a cursor program region, each qubit on the lower line has 
its own adjacent 4-state cursor system on the upper line 
that can be in one of the "active" states (J), ©, or Q. 
or in — , a new "inactive" state not used by V. Just as 
the program region qubits in V are always in classical 
states |0) or |1) for any simulation of a quantum circuit, 
so too will the qubits in the program regions of this ring 
always be in classical states, which we label schematically 
by the glyphs and 1 . However, qubits in the data 
region may be in any superposition of |0) and |1), and 
indeed even entangled with the state of all other data 
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qubits, so we represent an arbitrary state of a data qubit 
by the glyph ? . Table [I] is a legend for the glyphs we 
use to denote the cursor states, qubit states, and region 
boundaries in this geometry. 



B. Hamiltonian construction 

Like Feynman's Hamiltonian 0, we construct TJg on 
7^8 as a sum of terms with gi representing forward com- 
putation and g\ representing backwards computation: 
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FIG. 4: Glyphs for cursor states, qubits states, and region 
boundaries. (Color online.) 



Each operator gi is an interaction involving only the 8- 
state systems at sites i and (i + 1) mod K in the ring. A 
given g L is gcncrically responsible for two things. Firstly, 
it examines the current cursor state at site i, replaces it 
with — and appropriately sets the cursor at site i + 1 . 
This effectively passes the cursor down the ring and en- 
sures that only a single site is "active" (non- — ) . Sec- 
ondly, a gi is charged with performing any necessary gates 
on the qubits at sites i and i + 1, which depend on both 
the location i and the current state of the cursor. 

At most of the ring sites, gi merely swaps the state of 
the cursor. Namely, 
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where © represents an unspecified ®.® or © state 
that is preserved during the transition and ? represents 
either a (Q) or (V) that is preserved during the transition. 
In other words, Eq. Q is actually a sum of terms, one 
for each possible value of ? and ©. 

The only places where gi deviates from this behavior is 
near the boundaries between regions. This is because the 
circuit yW consists of SWAP gates everywhere except 
near these boundaries. We define how acts near these 
boundaries below: 
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• Swap Region Near the Swap-Data Boundary 

In order to simulate the action of V^ K \ which has 
a definite beginning and end, by Hg,, which has pe- 
riodic boundary conditions, we replace Eq. Q for 
the last two states in the swap program region with 
projections onto start and stop states: 
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FIG. 3: Architecture Layout (Color online) 



where the factors of 1/2 are to accommodate the 
fact that both gi and g\ are included in H$. 

The second term in ^ is a projection onto a stop 
state that is never reached by a valid programming 
of V^ K \ Specifically, "Behavior 4" in Sec. |h] pre- 
vents this from happening. (A cursor in the state 
© never sees a swap program qubit in the state 
©.) By augmenting the swap program with one 
extra qubit in the state 1 , the general swapping 
interaction ^ will evolve a state with a © atop 
© to the last two swap program sites discussed 
here and stop propagating further. Of course, be- 
cause a g\ term is in H$ for every gi term in Hg, 
it is possible for evolution to "undo" some of the 
computation and compute in reverse from the stop 
state. However, no additional forward computation 
will accrue after this stop state is reached. 

The first term in (|7| is also one that is never reached 
by a valid programming of V^ K \ Moreover, it is 
never reached by the extended "stop state" imple- 
mentation by Hg, just discussed. Indeed, it is pre- 
cisely the stop state projection term that prevents 
the start state from ever being reached. Hence no 
computation can be "undone" by g\ terms once this 
start state is reached (in reverse). 

The net effect of incorporating stop and start state 
interactions is that evolution by H% on a ring can be 
"unfurled" into an equivalent interaction by a re- 
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lated Hamiltonian on a line. This will be described 
in more detail in Sec. IIVI 

To simulate the action of on these two qubits 
for a valid programming of the swap program reg- 
ister, we need to add to gi of Eq. ^ the following 
terms: 



at the last two data sites with 
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+ 




®l\/® I 

©®IA®®I 

®i y ®ei 

®@l A®®l 



(8) 



Finally, although the following state only arises 
when V simulates one of the rather banal one-gate 
circuits U = (I <g) H)A(S) or U = A(S), we add to 
gi a projection onto a second stop state for com- 
pleteness' sake: 
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• Swap-Data Boundary 

To simulate the permutation of cursor states ® — > 
® — > O - ► ® enacted by when the state of 
the lower-most qubit in the swap program is a 0, 
we replace (J6j) at the swap-data boundary with 
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To simulate the movement of the cursor from the 
swap region to the data region when the last qubit 's 
state of the swap program is a 1 (and no cycling of 
the cursor's state is performed), we add to gi at 
this location one more term: 
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• Data Region Near the Data-Hadamard Boundary 

Because swaps the last two data qubits only 

when the cursor is in the state (§), and because 
y( K ) on ly applies a A(S) between these qubits 
when the cursor is in the state Q, we replace (JsJ) 
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where the bar over the qubits in the last two terms 
indicates that a A(S) gate has been applied between 
them. 

• Data-Hadamard Boundary 

Because no swapping occurs across the data- 
Hadamard boundary and because additionally a 
Hadamard gate is applied to the state of the last 
data qubit if the state of the upper- most Hadamard 
program qubit is a 1 and the cursor is in the 
state O we replace Eq. ^ at the data-Hadamard 
boundary with 
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where the bar over the qubit in the last term in- 
dicates that a Hadamard gate has been applied to 
it. 

• Hadamard- Swap Boundary 

The only subtlety about the Hadamard-swap 
boundary is that qubit states are not swapped 
across them in V^ K \ Hence Eq. ^ is replaced 
at this boundary by 
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Given this definition for our Hamiltonian Hg, our next 
task is to demonstrate how it may be programmed to 
simulate quantum circuits. In the next section, we show 
how this can be done when H% drives a continuous-time 
quantum walk. In the subsequent section, we show how 
this can be done when (a slightly modified version of) H% 
is the final Hamiltonian of an adiabatic algorithm. 
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IV. UNIVERSAL QUANTUM COMPUTATION 
VIA A CONTINUOUS-TIME QUANTUM WALK 
DRIVEN BY H 8 

The Hamiltonian Ha, defined in the previous section 
enables universal quantum computation by driving a 
continuous-time quantum walk in the following way. 
First choose a precision e and a quantum circuit W to 
simulate to within precision e. Then, construct a circuit 
U over the Kitaev gate basis that approximates W to 
within e. This can be done using standard techniques, 
e.g., by the Solovay-Kitaev algorithm (35]. Next, deter- 
mine how to program the circuit V defined in Sec. [TT]so 
that it simulates U exactly. Initialize the swap, data, and 
Hadamard regions of TLn®TL2 to the same states that one 
would initialize the swap, data, and Hadamard regions for 
V. However, because the swap region on the ID ring has 
one more state than in the V circuit, initialize the last 
swap qubit on the ID ring to be in the state |1). (Note 
that a valid programming obtained from V will necessar- 
ily also have the penultimate swap qubit also initialized 
to the state |1).) Finally, initialize the cursor line on the 
ID ring to be (3) everywhere, except over the last swap 
qubit, where it is in the state O We will denote the 
state of the ID ring so initialized by |-0o) - 

Only two gi terms in Ha, act nontrivially on \ipo). The 
first is the start state projector Q and the second is 
the transition term across the swap- data boundary (111. 



Thus evolution by H s can either keep the ring in the 
state \ipo) or advance it with some amplitude to a unique 
successor state \ipi). If the ring is in the state 
the Hamiltonian Ha, only couples it back to |^o) or to 
a unique successor state ^2) given by the general swap 
rule ([6]). This line of reasoning continues, with a unique 
successor state and predecessor state existing for \ip t ) in 
the ring until the final t = T := 0(nT). In this configu- 
ration, only the stop state projector and the predecessor 
state match. The Hamiltonian Ha, can therefore be re- 
stricted to a subspace of size T + 1 , on which it looks 
like 



ffr ) :=£l^>W>t-i| + K>t-i><^l (15) 

i=l 

+ l^)MI + W*r>(lM- ( 16 ) 

At this point, one can argue as Feynman originally did 
[15) that evolution by Ha, is the same as by a quantum 
walk on a line. Straightforward analysis demonstrates 
that the time at which this walk has a maximum ampli- 
tude for moving from \ip ) to is time T/2, at which 
the amplitude is 



-m 8 T/2 



T 



-1/3 



(17) 



©•■• 0001000010- -001 
©••• @®®i®®®®i©-- @®i 

FIG. 5: Start State: • ■ • represents an arbitrary number of 
qubit/qudit vertical pairs 



Given these insights, a continuous-time quantum walk 
driven by H s can be made to simulate V arbitrarily well 
in several different ways: 

• One could simply repeat the preparation and evo- 

2 /3 

lution 0(T ) times, availing upon the Chernoff 
bound that it is exponentially likely in T that one 
of the final measurements will find the system in 
the state representing the output of the com- 

putation. 

• One could replace the projection onto the stop state 

2 /3 

in H s with 0(T ) general swap transitions de- 
fined by Eq. ^ on C(T 2 ^ 3 ) additional qubits in the 
swap program region. By use of the Chernoff bound 
once again, it becomes exponentially likely in T 
that the cursor is in one of these new "dummy" lo- 

2 /3 

cations. Note that simply adding 0\T ) identity 
gates to the end of the circuit V being simulated is 
not trivial as it was in similar refs. |15l 137) because 
V simulates the identity by applying A(S') 4 . 

• One could replace both the start state and stop 
state projections in Ha, by additional "runway" and 

2 /3 

"landing pad" qubit chains of size 0(T ) that are 
swapped by the general swap rule ^ . As shown by 
Feynman [T5] (and argued in greater detail in |55] ) 
this will allow the computation to proceed ballisti- 
cally with arbitrarily high probability. 

• One could redefine Ha so that the tth term had 



A detailed analysis in terms of Bessel functions of why 
this is the case can be found in numerous places, for 
example in [36) . 



a coefficient y T(T — t) , enabling perfect fidelity 
state transfer to the final state of the computation 

HOH]. 

Running a quantum algorithm in this model is partic- 
ularly simple. One first initializes the state |^>o)> then 
waits a time T/2, and finally measures the state of the 
ring destructively using one of the methods described 
above. Because no dynamical controls are required dur- 
ing the course of the computation, "gate errors" appear 
as fabricational errors in Ha,. Rather than having to deal 
with such errors during runtime, which is computation- 
ally expensive and must be done in the circuit model, 
these errors may instead be dealt with during the fabri- 
cation of H$ in a much more controlled environment. De- 
coherence and other environmentally-induced noise pro- 
cesses will still be present, but at least they will not be 
conflated with dynamical control errors. 
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V. UNIVERSAL QUANTUM COMPUTATION 
VIA ADIABATIC EVOLUTION TO Hg 



Given that comparable physical models have appeared 
in proofs for adiabatic universality jST] H3J , we will 
demonstrate that one-dimensional Hamiltonians on 8- 
level quantum systems on a ring are universal under adi- 
abatic evolution. Our derivation will closely mirror these 
previous works and follows straightforwardly from results 
of Aharanov et al. [20] that have been further refined by 
Deift et al. [41]. We begin by quoting the adiabatic the- 
orem: 

The Adiabatic Theorem (adapted from |41j, 
quoted from [20j ) Let -ffinit and -ff nn ai be two Hamilto- 
nians acting on a quantum system and consider the time- 
dependent Hamiltonian H(s) = (1 — s)Hi nlt + s-ff nn ai- 
Assume that for all s, -ff(s) has a unique ground state. 
Then for any fixed S > 0, if 



t > n 



lltffi. 



nal 



init I 



\1+S 



1 sG[0,l 



{A^(H(s))} 



(18) 



then the final state of an adiabatic evolution accord- 
ing to H for time T (with an appropriate setting of 
global phase) is e-close in ^-norm to the ground state 
of -ff final- The matrix norm is the spectral norm ||ff|| = 
max w \\Hw\\/\\w\\. 

Loosely speaking, the adiabatic theorem states that if 
the variation between the Hamiltonians is slow enough, 
then the quantum state will be in the ground state of 
-fffinai at the end of evolution if it starts in the ground 
state of -f/initiai- Consequently, to demonstrate that our 
model is universal under adiabatic evolution, we must 
pick an .ff m it and -ffnnai whose ground states are appro- 
priately related to the initial and final states of compu- 
tation. 

First consider ffinit- We have previously detailed the 
configuration of a valid initial state, \tpo), which en- 
codes both the program of execution and the input data 
state. We define ff m it to be a sum over local projec- 
tors for which \ip ) is the zero-eigenvalue ground state. 
This differs from previous adiabatic proofs involving non- 
programmable architectures, which only needed to en- 
code the input data state. We will now construct -ffmit 
piece by piece. 

In a shorthand notation analagous to the one used in 



Sec. IV the first term is 




(19) 



where the projector is for the two locations to the left of 
the swap program/data region boundary. While this will 
ensure that the ground state is an initial state, it does 
not ensure that it is valid or that is executes the desired 
program (i.e. that it is \ipo))- However, we can pick 
out our desired initial state by adding nearest-neighbor 
projectors to this Hamiltonian in a chained fashion. The 



first such projector will be between the leftmost qubit 
in the above projector and its left neighbor. Without 
loss of generality, suppose that the program qubit for 
our desired program is a in this location; we would 
then add terms 



E 



©e v me 



E 





/me 




©n 1 


\@® 


ii 



(20) 
(21) 



where the first sum is over all qudit symbols, ensuring 
that anything over a (f) has a higher energy. Similarly, 
the second sum ensures that all configurations except the 
desired (3) over a (0) have a higher energy. Thus, we 
sum over all non-desired configurations, of which there 
are seven. This process is continued by sliding over one 
location and again summing over projectors where the 
right qudit-qubit pair is fixed to the desired initial con- 
figuration and the left qudit-qubit pair runs over all seven 
illegal or non-desired configurations. Moving around the 
entire ring step by step, we add similar projectors which 
tack on all undesired nearest-neighbor pairs. This en- 
sures that the ground state is precisely |^>o)j it is both 
valid and encodes only the program we want to execute. 

Now consider -ff nna i- As discussed in [201 [41] 02] , we 
seek a -ff nn ai whose ground state is the sum-over-histories 
state 



3=0 



(22) 



This can be accomplished simply by adding penalty 
terms to Eq. [15] which "prefer" the transition elements. 
That is, we define 



T-l 



-ffnnai := ^ \lpi) (ipi 



i = l 
T 



i=0 



(23) 

(24) 
(25) 



Expressed in the \ipt) basis, -ffmit is simply 



/o ... o\ 

1 ... 

V° ■•■ V 



(26) 
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as jV'o) is the unique groundstate. Similarly, -fffinai is 



o \ 



1 



- 







(27) 



As proved in Sec 3.1.2 of [2D] and simplified in [IT], the 
spectral gap of these Hamiltonians is at least 1/[2(T + 
l) 2 ], which is an inverse polynomial in T, the number 
of gates in the initial quantum circuit (recall T is poly- 
nomial in T). Thus, the evolution can be considered 
efficient, proving that Hamiltonians of 8-level quantum 
systems on a one-dimensional ring are universal for adi- 
abatic quantum computation. 



spatially-sparse fc-local Hamiltonian to a 2-local Hamil- 
tonian on a regular lattice. This reduction maintains the 
spatial setup of the initial Hamiltonian, so that for our 
ring geometry, the reduced setup would be a regular lat- 
tice on the surface of a cylinder. Moreover, the effective 
interaction is not only close to the initial Hamiltonian 
with respect to its eigenvalues, but the operator norm is 
close within a restricted subspace. Thus the dynamics are 
similarly close. In other words, given a desired e, a "re- 
duced" 2-local version of the architecture can be designed 
with dynamics which are e-close to that of the original 
architecture. While the actual mapping requires fine- 
tuning perturbation coupling parameters to ensure only 
polynomial overhead and growth, the Oliveira-Terhal ap- 
proach does provide a means for reducing the interaction 
degree of our architecture, though the resulting system 
is on the surface of a cylinder whose height is no longer 
fixed. 



VII. CONCLUSION 



VI. A TWO-DIMENSIONAL CYLINDER QUBIT 
HAMILTONIAN H 2 

Given the conceptual decomposition of the 8-level 
quantum system into a 4-level quantum system and a 
qubit, it is easy to imagine fully decomposing each 8- 
levcl quantum system into 3 qubits. Such an architecture 
would still be effectively one-dimensional, as the width of 
the ring is invariant. However, each nearest neighbor in- 
teraction in the 8-level system would become a 6-body 
interaction, which is a greater than Feynman's 4-body 
interactions. 

Alternatively, one could use gadget Hamiltonian the- 
ory, as presented in [19] , to reduce the degree of interac- 
tions between qubits. This procedure augments the orig- 
inal system S with additional mediator qubits. Given an 
initial Hamiltonian H acting on S, one constructs a new 
Hamiltonian H' which acts on the augmented system and 
whose action restricted to S is bounded as \\H' — H\ | < e 
for a desired e. For our purposes, we would seek to re- 
place the 6-body terms required for the transition rules 
with lower body terms. While in theory it is straightfor- 
ward to apply a 6-body gadget, each term in our original 
Hamiltonian will require 6 mediators. Generally, each 
6-body term will then require more than 6 mediators, 
resulting in an expanding geometry. We have been un- 
able to devise a scheme which maintains both spatially 
local interactions and an effective one-dimensional width 
under this replacement scheme. 

An alternative and more promising approach, pre- 
sented by Oliveira and Terhal [2Tj . details an "efficient" 
gadget reduction procedure to a 2-local Hamiltonian 
which is planar, but two dimensional. We defer to the 
paper for the details of such a reduction, but note that 
there is a prescription for applying gadgets to reduce a 



We presented a family of time-independent Hamilto- 
nians that can enable universal quantum computation 
either by driving a continuous-time quantum walk or 
by terminating an adiabatic algorithm. When used to 
drive a continuous-time quantum walk, quantum compu- 
tation consists of 1) preparing an input that describes a 
quantum circuit to be executed and the quantum data 
onto which it is to be applied, 2) waiting the appropriate 
amount of time, and 3) measuring the output. The sim- 
plicity of the operation of such a machine is appealing, 
but Feynman's original proposal for realizing it [T5] still 
remains far out of technological reach. Our work demon- 
strates that it suffices for such a machine to use only 
nearest-neighbor interactions between 8-level systems on 
a ID ring, which may be more technologically feasible. 
It also demonstrates that simulating the dynamics of ID 
time-independent Hamiltonians on 8-level systems is a 
BQP-complete problem, even though it is known how to 
simulate ID spin systems in a way that scales polynomi- 
ally with the number of spins [43] . 

When used to terminate an adiabatic algorithm, our 
Hamiltonian achieves universality via nearest-neighbor 
interactions between 8-level systems on a ID ring, one 
level fewer than in a recent proof of universality of 
the adiabatic algorithm by 9-level nearest-neighbor ID 
Hamiltonians [2"3"] . 

Finally, using gadget perturbation theory [21] , our 
Hamiltonian can be made spatially local using only 
qubits rather than 8-level systems, but at the expense of 
requiring the qubits to lie on the two-dimensional geome- 
try of a cylinder rather than a one-dimensional geometry 
of a ring. 

A remaining challenge is to address error correction 
and fault tolerance for adiabatic and quantum walk mod- 
els driven by the Hamiltonian we constructed [44] . In 
particular, because the quantum walk describing quan- 
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turn computation is effectively a continuous-time quan- 
tum walk on a line, an imperfect implementation may be 
subject to Anderson localization which could exponen- 
tially suppress propagation along this line. One promis- 
ing feature of the quantum walk model is that there are 
no dynamical controls during the operation of the quan- 
tum computer, so that control errors are all fabricational. 
This allows such errors to be handled at "compile time" 
rather than at "run time," which could be much easier. 

Another remaining challenge is to explore whether 
finding the ground state of our Hamiltonian is a QMA 
complete problem. We conjecture that it is, as several 
similar Hamiltonians used to prove adiabatic quantum 
computing universality are j2Qj [23] . However, the "clock" 
and its update rule as implemented by our Hamiltonian 
are not localized to a particular point in space, as com- 
putation winds around the ring many times during the 
course of a computation. Hence, previous proofs do not 
translate directly and creative ideas are required. Since 
our architecture is also programmable, the need to en- 
code the program in the initial state further complicates 
the QMA question. 



Addendum: As we were finishing this paper, we 
became aware of related work by Nagaj and Wocjan 
that demonstrates that translationally invariant quan- 
tum walks in one dimension ("continuous-time quantum 
cellular automata" ) are also universal for quantum com- 
putation, albeit using ten-dimensional rather than eight- 
dimensional systems |32j . A final remaining challenge 
we state is to explore whether our model can be made 
translationally invariant or whether the dimension of the 
Nagaj- Wocjan model can be reduced to eight as ours is. 
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APPENDIX A: EXAMPLE COMPUTATION 



The following illustrates computing the circuit in Figure [6] 




FIG. 6: Equivalent Circuit 



11 



Start state 

Swap data 

A(5) on ? # 

-ff since control bit is (T) 

Cycle Hadamard 



Ol 901 06 -I 



2. 
3. 
4. 
5. 



000001 0001 00001 
@@@®0 1 1 ?1®@®(g 

©00001 0O01 00001 

@@@®o 1 1 1 m 1 



00000 



N ®®@(§ 



! 00IO000I 

' l®@(oXo 



big 3\Qoe 

©©©00 o®@@l 



! 9I00O0I 

H@@®@l 



Cycle swap 9. 

10. 
11. 
12. 

Become (J) across boundary 13. 
Cycle data 14. 

15. 



310 

@@@®®l 

O0000IG 

@@@®®l 



0©0£ 010 301 ©0£ 

@@@®C o©@®l 



30G I 0100001 
H@@@®l 



30GI 0100001 

H@@@®l 



©000OI0 3i eeeei 

@@©®(o| I I •! '@@®l 



000001® 

@@®®@l 



000001 0®01 0000 

@@®®(o I I ? •! o @(Q)® 

00000100®! 0000 
@@®®(o| ?•? |o @@®l 



Cycle Hadamard 



17. 



0000010001 ®00©l 
@@®®(0 | I • I I (Q)@®l 



18- 300001 
@@®®@l 



0I0®00 
1 1 o)@@® 



19. 



3IG30I00®G 
@@®®(o I ? # ? I o<§Xg)Ci 



Cycle swap 



21. 
22. 
23. 
24. 

Become ® across boundary 25. 



310001000®! 
@@®®(o| ?®9!@©®(g 
©300010001 00001 

1 1 0| ? # ?| 1 0| 



)®0©0K 

1101 



01000© 

? I oj©(D@ 



00®00I000IG 
@®@®(Q| ?fy. 



®00®( 
@®® o o 



100001 

? | 0)@®@| 



3®IG 

@®®:o Q| Q'@®@l 



12 



Swap data 26. 9|®G H 

@® ® @@l @@®@l 

Do not swap data 27. 

@®®@(Q|#? ? | oxo)®@l 

28. ©0©00|00®|0000| 

@(D(D(oXo|#? ? | o)@®@l 

Cycle Hadamard 29. 000001 0001 ©0001 

@®(D(0)0|+? ?|0)@®@l 

30. 00000|000|0®00| 
@®®(§)0|»? ?|0)@®@l 

31- 000001000100001 
@®®@'0|»? ? I 1 I 

32. 00000|000|000®| 

@®®(0)0|+?»l@®@(Q)| 

Cycles 33. ®- 9|®( 

@®®(0>0|#? »l@©@@l 

34. 0®000|000|0000| 
®@®@'0|#? ?|o)®@(o)l 

35. ©0©00|000|0000| 
®®@@'0|#? ?|0)®@@l 

36. 000®0|000|0000| 
®®@@(0|+?»l@®@@l 

Become O across boundary 37. 

®®@@@l» | ? I 0)®@@l 

Swap data 38. 9|QG 

®®@@0|»?»l@®@(0)| 
A(S)on#? 39. 3000|0Q0|00 

®®@(0>0|?»?|0)®@(0)| 

No H since control bit is 40. 9|£ 

®®@@(0 I I #T I ®®®! 
Cycle Control Program 41. j 9|©G 

® ® ® ® ®l ® ® ®@ I 

Become Q across boundary 42. 9|£ 9|®OG 

®®@@(o i ? m i|®@@@| 

43. 00000I000I00O0I 
®®@®(0| ?#»|®@ 0| 

44. ©0000I000I000OI 
®®@®:o| V»l®@®@l 

Cycle swap 45. Q 9©0©| ©0©| 0©0©| 

®® ® ® ®l ® ® ®® I 

46. 0O00©I000I000©I 
® ® ® ® ®l ® ® ®® I 

47. 00O00I000I000G 
® ® ® ® ®| •••! ® ® ®® I 

Stop state 48. 3O0IG 

® ® ® ® ®l ®@®@l 



[1] D. Deutsch, Quantum computational networks, Proc. [2] E. Bernstein and U. Vazirani, Quantum complexity the- 
Roy. Soc. London A 425, 73 (1989). 



13 



ory, in Proceedings of the 25th Annual ACM Sym- 
posium on the Theory of Computing, ACM (ACM 
Press, New York, 16-18 May 1993, San Diego, 
CA, USA, 1993), pp. 11- 20, ISBN 0-89791-591-7, 
|doi : 10 . 1145/167088 . 167097] 
[3] A. C.-C. Yao, Quantum circuit complexity, in Proceed- 
ings of the 34th Annual Symposium on Foundations 
of Computer Science, IEEE (IEEE Press, New York, 
3-5 Nov. 1993, Palo Alto, CA , 1993), pp. 352-361, 
|doi : 10 . 1 109/SFCS ■ 1993 ■ 366"852l 
[4] A. Barenco, A universal two-bit gate for quantum com- 
putation, Proc. Roy. Soc. London A 449, 679 (1995), 
arXiv : quant-ph/9505016 
[5] D. P. DiVincenzo, Two-bit gates are univer- 
sal for quantum computation, Phys. Rev. A 51 
1015 (1995), |doi : 10 . 1103/PhysRevA . 51 . 1015 



|arXivTc ond-mat/9407022r 
[6] S. Lloyd, Almost any quantum logic gate is 
universal, Phys. Rev. Lett. 75, 346 (1995), 
|doi : To". 1 103/PhysRev Lett . 75 . 346] 
[7] P. W. Shor, Fault-tolerant quantum computation, in Pro- 
ceedings of the 37th Annual Symposium on Foundations 
of Computer Science, edited by R. S. Sipple, IEEE (IEEE 
Press, Los Alamitos, CA, 14-16 Oct. 1996, Burling- 
ton, VT, USA, 1996), pp . 56-65, ISBN 0-8186-7594-2, 
|arXiv : quant -ph/960501 1] 
[8] A. Y. Kitaev, Fault-tolerant quantum computation by 

anyons (1997), [arXiv: quant -ph/ 9707021 
[9] D. Aharonov and M. Ben-Or, Fault tolerant quantum 
computation with constant error rate (1999), See also 
condensed version [45]., arXiv : quant-ph/9906129 

[10] P. O. Boykin, T. Mor, M. Pulver, V. Roy- 
chowdhury, and F. Vatan, On universal and 
fault-tolerant quantum computing, pp. 486- 
494 (1999), |doi:0.1109/SFFCS. 1999.814621 

|arXiv : q uant-ph/9906054] 

[11] P. O. Boykin, T. Mor, M. Pulver, V. Roychowdhury, and 
F. Vatan, A new universal and fault tolerant quantum 
basts, Inf. Proc. Lett. 75, 101 (2000). 

[12] Y. Shi, Both Toff oh and controlled-NOT need little help to 
do universal quantum computation, Quant. Info. Comp. 
3, 84 (2003), [arXiv : quant-ph/0205115l 

[13] R. Solovay, Lie groups and quantum circuits (2000), 
MSRI presentation at http://www.msri.org/ 



publica-tions/ln/msri/2000/qcomputing/ solovay/ 



3zr 



[14] A. Y. Kitaev, Quantum computations: algorithms and 
error correction, Russian Math. Surveys 52, 1191 (1997). 

[15] R. P. Feynman, Quantum mechanical computers, Found. 
Phys. 16, 507 (1986). 

[16] E. Farhi and S. Gutmann, Quantum compu- 
tation and decision trees, Phys. Rev. A 58, 

915 (1998), |doi : 10 . 11 03/PhysRe vA . 58 ■ 915| 

arXiv : quant-ph/9706062 

[17] R. Landauer, The physical nature of information, Phys. 
Lett. A 217, 188 (1996). 

[18] A. Yu. Kitaev, A. Shen, and M. N. Vyalyi, Classical 
and Quantum Computation, vol. 47 of Graduate Studies 
in Mathematics (American Mathematical Society, Prov- 
idence, RI, 2002), ISBN 0-821-82161-X, Translated from 
the Russian edition by Lester J. Senechal. 

[19] J. Kempe, A. Kitaev, and O. Regev, The complexity of 
the local Hamiltonian problem, SIAM J. Comput. 35, 
1070 (2006), [arXiv: quant -ph/0406180l 



[20] D. Aharonov, W. van Dam, J. Kempe, Z. Landau, 
S. Lloyd, and O. Regev, Adiabatic quantum com- 
putation is equivalent to standard quantum computa- 
tion, in [37], pp. 42-51 , |doi : 10 . 1 109/FOCS ■ 2004 . 8| 
arXiv : quant-ph/0405098] 

[21] R. Oliveira and B. M. Terhal, The complexity of quantum 
spin systems on a two-dimensional square lattice (2005), 
|arXiv : q uant-ph/0504050] 

[22] K. G. H. Vollbrecht and J. I. Cirac, Quantum simulators, 
continuous-time automata, and translatioanlly invariant 
systems (2007) , [arXiv : 0704 . 3432] 

[23] D. Aharonov, D. Gottesman, S. Irani, and J. Kempe, 
The power of quantum systems on a line (2007), 
larXiv: 0705 .40771 

[24] R. A. Horn and C. R. Johnson, Matrix analysis (Cam- 
bridge University Press, Cambridge, 1985), ISBN 0-521- 
30586-1. 

[25] N. Margolus, Parallel quantum computation, in Com- 
plexity, entropy, and the physics of information, edited 
by W. H. Zurek (Addison- Wesley, Redwood City, CA, 
1990), vol. VII of SFI Studies in the Sciences of Com- 
plexity, pp. 273-287, ISBN 0-201-51509-1, UR L |http:| 
I //kh . bu . edu/qcl/m argolunl9901clf 6f 6d . pdf 

[26] D. ~K. Meyer, From quantum cellular automata to 
quantum lattice gases, J. Stat. Phys 85, 551 (1996), 
|doi : 10 . 1007/BF02199356| [arXiv : quant-ph/9604003| 

[27] W. van Dam, A universal quantum cellular automaton, 
in Proceedings of the Fourth Workshop on the Physics 
of Computation, edited by T. Toffoli, M. Biafore, and 
J. Leao (New England Complex Systems Institute, Cam- 
bridge, MA, 22-24 Nov. 1996, Boston, MA, USA, 1996), 
pp. 323-331, URL http://www.cs.ucsb.edu/~vandam/ 
research/uqca.pdf] 

[28] J. Watrous, On one- dimensional quantum cellular au- 
tomata, in Proceedings of the 36th Annual Symposium 
on Foundations of Computer Science, IEEE (IEEE 
Press, Los Alamitos, CA, 23-25 Oct. 1995, Milwau- 
kee, WI, USA, 1995), pp. 528-53 7, ISBN 0-8186-7183-1, 
|doi : 10 . 1 109/SFCS ■ 1 995 . 492583] 

[29] D. J. Shepherd, T. Franz, and R. F. Werner, 
Universally programmable quantum cellu- 
lar automaton, Phys. Rev. Lett. 97, 020502 

(2006), |doi : 10 . 1 103/PhysRevLett . 97 ■ 020502 

|arXiv:quant- ph/0512058] 

[30] D. Janzing and P. Wocjan, Ergodic quan- 
tum computing, Quant Inf. Proc. 4, 129 

(2005), |doi : 10 . 1007/slll28-005-4482-9| 

arXiv : quant-ph/0406235 

[31] D. Janzing, Spin-1/2 particles moving on a 2d 
lattice with nearest-neighbor interactions can realize 
an autonomous quantum computer, Phys. Rev. A 
p. 012307 (2007), |doi:10.1103/PhysRevA.75.012307| 
|arXiv : q uant-ph/0506270] 

[32] D. Nagaj and P. Wocjan, Hamiltonian quantum cellular 
automata in ID (2008) , arXiv : 0802 . 0886 

[33] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000), ISBN 0-521-63235-8 (Hardback), 0- 
521-63503-9 (Paperback). 

[34] A. M. Turing, On computable numbers, with an 
application ot the Entscheidungsproblem, Proc. 
Lond. Math. Soc. 42, 230 (1 937), errata in [46]., 
|doi : 10 . 1 112/plms/s2-42 . 1 . 230| 

[35] C. M. Dawson and M. A. Nielsen, The Solovay- 



14 

Kitaev algorithm, Quant. Inf. Comp. 6, 81 (2006), 
arXiv : quant-ph/0505030 
[36] A. M. Childs, R. Clove. E. Deotto, E. Farhi, S. Gut- 
mann, and D. A. Spielman, Exponential algorithmic 
speedup by a quantum walk, in Proceedings of the 
Thirty-Fifth Annual ACM Symposium on the The- 
ory of Computing, edited by L. L. Larmore and 
M. X.' Goemans, ACM (ACM Press, New York, 9- 
11 Jun. 2003, San Diego, CA, USA, 2003), pp. 59- 



68, ISBN 1-58113-674-9, 


doi : 10 . 1145/780542 . 780552 


arXiv:quant-ph/0209131 





[37] D. Azada, ed., Proceedings of the 45th IEEE Symposium 
on Foundations of Computer Science, IEEE (IEEE Press, 
Los Alamitos, CA, Oct. 17-19, Rome, IT, 2004), ISBN 
0-7695-2228-9. 

[38] T. J. Osborne and N. Linden, The propagation of quan- 
tum information through a spin system, Phys. Rev. A 
69, 052315 (2004), |doi: 10. 1103/PhysRevA.69.052315| 
arXiv : quant -ph/0312 141| U RL |http : //link . aps ■ org/| 
abstract/PRA/v69/e052315 

[39] M. Christandl, N. Datta, A. Ekert, and A. J. 
Landahl, Perfect state transfer in quantum 
spin networks, Phys. Rev. Lett. 92, 187902 

(2004) |doi:10. 1 103/PhysRevL ett . 92 . 187902 

|arXiv : q uant-ph/030913l] 

[40] M. Christandl, N. Datta, T. C. Dorlas, A. Ekert, 
A. Kay, and A. J. Landahl, Perfect transfer of arbi- 
trary states in quantum spin networks, Phys. Rev. A 
71, 03 2312 (20 05), |doi : 10 . 1103/PhysRevA . 71~ 032312 
|arXiv : quant -ph/041 1020] 

[41] P. Deift, M. B. Ruskai, and W. Spitzer, Improved 
gap estimates for simulating quantum circuits by adi- 
abatic evolution, Quant. Inf. Proc. 6, 121 (2007), 
|arXiv : quant - ph/0605 156] 

[42] S. Irani, The complexity of quantum systems on a one- 
dimensional chain (2007), arXiv: 0705. 4067 

[43] T. J. Osborne, The dynamics of Id quantum spin sys- 
tems can be approximated efficiently, Phys. Rev. Lett. 97, 
157202 (2006) , |doi : 10 . 1 103 /PhysRev Lett . 97 ■ 157202] 
|arXiv : quant -ph/050803i"1 

[44] D. A. Lidar, Towards fault tolerant adiabatic quantum 
computation (2007) , [arXiv : 0707 . 002l| 

[45] D. Aharonov and M. Ben-Or, Fault tolerant quan- 
tum computation with constant error, in Proceedings 
of the Twenty-Ninth Annual ACM Symposium on the 
Theory of Computing, edited by F. T. Leighton and 
P. Shor, ACM (ACM Press, New York, El Paso, TX, 
USA, 1997), pp. 176-188, ISBN 0-89791-888-6, See 
also extended version [SJ., doi : 10 . 1145/258533 . 258579 
|arXiv : quant -ph/9 61 1025 

[46] A. M. Turing, On computable numbers, with an appli- 
cation ot the Entscheidungsproblem. a correction., Proc. 
Lond. Math. Soc. 43, 544 (1937). 



